BANACH CENTER PUBLICATIONS, VOLUME ** 
INSTITUTE OF MATHEMATICS 
POLISH ACADEMY OF SCIENCES 
WARSZAWA 200* 



THE WEYL ALGEBRA, SPHERICAL HARMONICS, 
AND HAHN POLYNOMIALS 

EWA GNATOWSKA 
Department of Mathematical Methods of Physics, Faculty of Physics, University of Warsaw 

Hoia 74, 00-682 Warsaw, Poland 
E-mail: gnatowsk@fuw. edu.pl 

ALEKSANDER STRASBURGER 
Institute of Mathematics, University of Bialystok 
Akademicka 2, 15-267 Bialystok 
and 

Department of Mathematical Methods of Physics, Faculty of Physics, University of Warsaw 

Hoia 74, 00-682 Warsaw, Poland 
E-mail: alekstra@math.uwb.edu.pl 



Abstract. In this article we apply the duality technique of R. Howe to study the structure 
of the Weyl algebra. We introduce a one-parameter family of "ordering maps", where by an 
ordering map we understand a vector space isomorphism of the polynomial algebra on R 2d 
with the Weyl algebra generated by creation and annihilation operators oi, . . . , a^, a^~, . . . , . 
Corresponding to these orderings, we construct a one-parameter family of sl2 actions on the Weyl 
algebra, what enables us to define and study certain subspaces of the Weyl algebra — the space 
of Weyl spherical harmonics and the space of "radial polynomials" . For the latter we generalize 
results of Luck and Biedenharn, Bender et al., and Koornwinder describing the radial elements 
in terms of continuous Hahn polynomials of the number operator. 

1 Introduction. This article reports on an investigation (still in progress) of the 
Weyl algebra with the aid of R. Howe's duality technique. More specifically, we attempt 
to transfer the information obtained by considering the action on the (commutative) 
polynomial algebra on R 2rf of a triple of classical operators (consisting of the Laplacian, 
Euler operator and multiplication by the square of the radius and giving rise to an 
SI2 action on this algebra) to the case of the (noncommutative) Weyl algebra with 2d 
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generators — see the beginning of Section 3 for the definition of the Weyl algebra. On 
this route one is inevitably confronted with the renowned "ordering problem of quantum 
mechanics" , c.f. e.g. [[l7| |], || . We deal with this problem in a wider context by introducing 
a one-parameter family of "orderings" , by which we understand a family of vector space 
isomorphisms of the polynomial algebra with the Weyl algebra. Corresponding to these 
orderings, we construct a one-parameter family of 0I2 representations acting on the Weyl 
algebra, whose action we employ to decompose the algebra. 

Two subspaces of the Weyl algebra stand out in our analysis — the subspace of 'Weyl 
spherical harmonics" and the (in a certain sense complementary) subspace of "radial" 
polynomials, which are polynomials of a single element, the so-called number operator. 
In this paper we concentrate on the latter one, emphasizing the connection with the 
problems of the special function theory. This is the subject of the second part of the 
paper, where we describe properties of these radial polynomials in various orderings. In 
particular we identify the elements of a natural basis of the space with Hahn polynomials 
of the shifted number operator. This is a generalization of results obtained earlier by 
Lohe, Biedenharn and Louck pl[ , Bender, Mead, and Pinsky [^J, and Koornwinder 
among others. 



2 Preliminaries on the polynomial algebra in even dimension. Consider the 
real cartesian space of even dimension R 2d , d > 1, with coordinates denoted by x\, . . . , 
x d, £1, • ■ ■ j £,d- Let P = Pc(M. 2d ) denote the polynomial algebra with complex coefficients 
over K 2d and P k = P^(IR 2d ) its subspace consisting of homogeneous polynomials of degree 
k. Setting Zj = Xj + i^j and Zj = xj — i^j for j = 1, . . . , d, where here and everywhere in 
the sequel i = V— T, we can regard P as the polynomial algebra with respect to coordinate 
functions z\ 1 . . . , Zd and their conjugates z\, . . . , z~d- Thus, using the usual multi-index 
notation we shall write its elements in the form 

(2.1) p(z,z) = Y,P^ za ^> 

a,/3 

where z — (zi, . . . , Zd), z — (fi, . . . , z~d) £ C rf and p a ^ € C. In particular, 

d d d 

is the square-of-the-radius function (euclidean length squared) on R 2d . 
Similarly, setting as usual 

d _ 1 / d . d \ d _ 1 / d .d_ 

dzj 2 \ dxj d£j J ' dzj 2 \ dxj d^j 
one can express in the complex form the euclidean Laplace operator (Laplacian) 

d Q2 d Q2 d q2 

A = ^arf + ^M = 4 ^ dzjdzj ' 

J = l 3 J = l ^3 = \ J J 

and the Euler operator 

' '' d d c) d 



j=i j j=i ^ j =1 j j =1 
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Along with the usual action of the orthogonal group 0(2d) on P, we shall be concerned 
with the action of the unitary group U(e2) on P which is obtained by extending the action 
on generators given by 

z^gz, z^gz, g £ U(«f), 

where z — [z\, . . . , Zd) and z — {z\, . . . , Z4) are regarded as column vectors. Recall that 
the contragredient matrix (<7*) _1 of any g € U(d) satisfies (<7*) _1 = g, the bar denoting 
complex conjugation (taken entry- wise). These two actions are related by the natural 
embedding of U(d) into 0(2d). 

Modifying slightly the above differential operators and supplementing them by the 
multiplication by the square of the radius one obtains a triple of endomorphisms inducing 
an action of the Lie algebra sfe on P. In fact, it can be easily checked that by setting 

(2.2) Rp = r 2 Pl Lp=^Ap, Ep=(E + d)p, for p e P, 

(here E + d is the symmetrized Euler operator) we obtain commutation rules for the 
standard generators of the Lie algebra SI2, i.e. 

(2.3) [R, L] = -E [E, R] = 2R [E, L] = -2L. 

It can be verified that this latter action commutes with the action of 0(2d). A remark- 
able although elementary application of Howe's duality principle based on using the above 
actions is the analysis of the structure of the algebra P (cf. [||, p. 118), which gives a large 
part of the classical theory of spherical harmonics. In particular one arrives this way at 
the decomposition P = J ®H, where J C P is the subalgebra of 0(2d) invariant polyno- 
mials, known to be the subalgebra generated by r 2 , and H = H(R 2d ) = {p E P \ Ap = 0} 
is the subspace of harmonic polynomials in P. The isomorphism is obtained via multi- 
plication. Consequently each homogeneous polynomial p £ P can be decomposed into a 
sum of products of homogeneous harmonic polynomials with powers of the radius-square; 

[3 dc SP] 

P(z,z) = ^2 r^hji 2 ,*), 
3=0 

where hj are homogeneous harmonic polynomials of degree degp — 2j . (See e.g. [^j for a 
derivation of this decomposition along classical lines). 

Furthermore, the following facts can be obtained within this approach. 

Let H k C H denote the subspace of homogeneous harmonic polynomials of degree k. 

a) For each k £ N the subspace H k is irreducible under the natural action of 0(2g?) 
and for any j ^ k the actions on £P and H k are not equivalent. Therefore the direct sum 
decomposition H — 0°^ o i? J is multiplicity free. 

b) For each to 6 N the space P m decomposes into a direct sum of 0(2d) irreducible 
subspaces as 

[im] 

(2.4) P rn = © r 2j ir"- 2j . 

c) Let P(R 2d ; H k ) be the sum of all 0(2d) irreducible subspaces of P equivalent to 
H k ; it is called isotypic subspace of type H k . Then the multiplication 

(2.5) J <g> H k -> P{R 2d ;H k ); p(r 2 ) ® h(z, z) h-> p(r 2 ) ■ h(z, z) 
induces an equivariant isomorphism of P(M. 2d ; H k ) with J® H k . 
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As the actions of and 0(2d) commute with each other, it makes sense to consider 
their joint action — perhaps more properly one should speak of the action of the direct 
product of the Lie algebras and o(2d) on P. It can be shown that the space P(M. 2d ; H k ) 
is irreducible under this joint action. By transfering the action of to the algebra of 
invariants J = P[r 2 ] by means of the multiplication map (|]^) one obtains the action 
which is (algebraically) equivalent to that of one of the lowest weight representations of 
5I2 — the lowest weight (with respect to E) being equal to (k + d). Actually, the lowest 
weight vectors in this space can be identified as certain orthogonal polynomials. 

Lastly, it is known, that under restriction to U(c?) C 0(2d) the irreducible 0(2c?)-space 
H k splits into the direct sum of U(ci)-irreducible subspaces, 

(2.6) H k = H {n ^ 

m-\-n—k 

where H^ n ' m ' denotes the space of homogeneous polynomials of bi-degree (n, to). We 
recall that p(z, z) £ fl"( n > m ) <^=> p(Xz, \z) = X n \ m p(z, z). Consequently, the isotypic 
subspaces of type if( n ' m ) in P are of the form J • H^- n ' m \ and 

P(R 2d ;H k ) = J-H {n ' m \ 

m+n—k 



For the explicit procedure of obtaining the decomposition (2.6), in particular for the 



formulae expressing projections onto subspaces of the given bi-degree, we refer to 



3 Construction of an sfe action on the Weyl algebra. In this section we shall 
describe a method of transfering the results concerning the polynomial algebra described 
at the end of the preceeding section to the case of the Weyl algebra. 

Usually, by the Weyl algebra with 2c? generators W = W2d one means the associa- 
tive algebra with unit I generated by 2d elements eii, . . . , ey, , . . . , al" subject to the 
relations (called Canonical Commutation Relations — abbreviated CCR in the sequel) 

(3.1) ajO^ - a^a 3 = 5 kj I, a,a fe - a k aj = 0, a+ - a^a+ = 0, 1 < j,k < d. 

We shall employ the following standard notation from the theory of associative al- 
gebras: for arbitrary x, y € W we shall denote by [x, y] — xy — yx the commutator of 
x and y and by {x, y} — xy + yx their anticommutator. For x € W the adjoint map 
ad(s) : W — * W is defined by setting ad(x) w — [x, w] for w £ W — as well-known, 
ad(x) is a derivation of VV. 

Let us recall that for many questions it is more convenient to realize the abstractly 
defined Weyl algebra as the algebra of partial differential operators with polynomial 
coefficients, considered initially as operators acting as on the Schwartz space <S(IR d ), 
and regarded, when need arrises, as densely defined operators in L 2 (R d ,dA(x)) — the 
space of square integrable functions on R d with respect to the Lebesgue measure. This 
identification is obtained by the standard assignment 

Note that a* = a^~, where the adjoint * operation is understood in the sense of the formal 
adjoint of a differential operator. In quantum physics the operators aj, are known as 
annihilation and creation operators, respectively. For the present work however, it is the 
algebraic structure of the Weyl algebra W that matters. 
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Given a 6 VV we denote by L a the map VV — *■ W obtained by the left multiplication 
with a and by the map corresponding to the left multiplication by a + . The symbols 
R a , R£ will have analogous meanings with respect to the right multiplication. By tak- 
ing convex combinations of left and right multiplication operators we shall be able to 
construct a family of endomorphisms of VV which will serve as a replacement for the op- 
erators of multiplication by coordinates and their complex conjugates in the polynomial 
algebra. 

Definition 1. For a given a <G VV and q 6 [0, 1] we define (a) endomorphisms M a and 
M+ of VV by the formulae 



(3.2) M a = (1 - q)L a + qR a = L a -q ad(o) 

(3.3) M+ = qL+ + (1 - q)R+ = L+ - (1 - q) ad(a+) . 

Let us explicitly point out to the meaning of this definition for the special values 
q = 0, |, 1 of the parameter. For (7 = 0, M a is the left multiplication by a, while M+ is 
the right multiplication by a + . For q = 1 these roles are reversed, while for q = \ each is 
equal to 1/2 of the respective anticommutator. 

Lemma 1. Let {ai, . . . , a^, a^, . . . , aj~} be a set of generators of the Weyl algebra 
and put Mj = M aj , M+ = M+ , j = 1, . . . , d. Then {Mi, M d , , . . . , Mf} is a 
commutative family of endomorphisms ofW. 

Proof. In view of CCR it is clear that the Mj's (M^'s respectively) commute between 
themselves, and that [Mj, M£] = for j ^ k. This last relation is also true for j = k. In 
fact 

[M„M+] = (1 - q)q([L aj , L+] + [R aj , fl+ ]) = (1 - g)?(I - I) = 0.D 
Thus employing the multi-index notation for the family {Mi, . . . , Md, M±, . . . , Mj"}, 
i.e. writing M a M +f} = Mf 1 . . . M* d M+ 01 . . . M+ Pd for a,/5eZf,we have by virtue of 
the lemma 

(3.4) M a M +f3 ■ M a 'M +fJ ' = M {a+a 'Hl + ^ +fi '\ 



Lemma 2. The map O q : P(R 2d ) -> W 2d de/med 6y 

(3.5) O q {z a z^) = M a M +f3 I 
is a vector space isomorphism such that 

(3.6) O q (z jP ) = Mj-O,^), 0,(z#) = M+O q {p), 

(3.7) G 9 f = ad(a+) 0,(p), 9 ( %L) = adfo) C,(p). 



Proof. The first pair of relations is obvious by the definition of O q and Lemma 1. The 
second pair, concerning derivations, follows by routine calculation. □ 

Definition 2. The map O q : P(R 2d ) —> W2d defined for q £ [0,1] by the equation 
(3.5) will be called the g-ordering map. 



1 The condition q G [0, 1] turns out to be nonessential and will be removed in further 
publications. 
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Since the subalgebras of W generated by pairs {a/, a z + }, {at, a^} with k ^ I, commute 
with each other, the elements ([T^) can be written in the factorized form, 

d 

O q (z a &) = l[B%(a l) a+), 

1=1 

where for nonnegative integers k, j we have set 

k 

(3.8) 



In the special cases mentioned above the formula (3.8) specialises to the following 
well-known constructions. 

a) For q = it gives the so-called normal (Wick) ordering 

O (z a zP)=a a (a+f, 

with all a,j (annihilation operators) on the left and all a~j (creation operators) on the 
right. 

b) For q = 1 the roles of dj and are reversed and 0\ gives the anti-normal 
(anti-Wick) ordering 

Oi(z a zP) = a +(} a a . 

c) The case q = h corresponds to the so-called Weyl, or symmetric, ordering. 



All these assertions are evident on the basis of (3.2). Further, let us observe that our 
elements Bj are related to the basis elements introduced for the case of one degree 
of freedom (d = 1) by Louck and Biedenharn in ||. In fact, Bj is proportional to 



p' J S(a ) a + ) = 2 L ^. 



'k\ ^ a l (a+ya k - 1 



Since the q-ordering map O q is for each value of q an isomorphism of vector spaces, there 
exist cndomorphisms of the Weyl algebra, which we shall denote lZ q , C q , and £ q (for 
simplicity we suppress the index q writing 1Z = 7i q etc., when no confusion can arise), 
such that the following equalities hold 

(3.9) O q R = K q O q , O q L = C q O qi O q E = £ q O q . 

This can be conveniently summarised in the form of commutative diagrams 

P c (R 2d ) "- Ll< P c (R 2d ) 



(3.10) 



w 2d 



T^q j C-q ; £q 



w 2d 



^From (2.2) and (p. 6 p. 7) we can deduce explicit expressions for lZ q , C q and £ q . 

Lemma 3. The endomorphisms lZ q , C qi £ q of the Weyl algebra W 2 d are given by the 
following formulae: 

d d d 

(3.11) lZ q w = (1 — q) 2 ajwaf + g(l - q) \2( wa f a j + a j a j w ) + I 2 afwdj, 
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(3.12) C q w — — ad(aj) ad(a^ 

(3.13) f,w = — / — q)(aj[af,w] — [a 3 -, w]af) + g([aj", tu] a j ~ a tl a ji w ])} + ^ 

3=1 

Note that C q does not depend on q and moreover, since ad(aj) ad(a^") = ad(a^~) ad(aj), 
we also have C q = — Ylj=i ad(a^) ad(<Zj). 

For the case of the Weyl (symmetric) ordering, i.e. for<j= i, we shall employ the no- 
tation 1Z S , C s , £ s for these operators. Note that they are given by fairly simple expressions 
(recall {a,b} stands for the anticommutator of a, b £ W): 

1 d 

j=i 

d 

£ s w = -^2[aj,[aj,w\], 
j=i 

d 

£ s w = ^(ajwoj" — a~j~waj). 



Applying successively the defining relation (|3.9| ) to the commutation relations (2.3), we see 
that the endomorphisms lZ q , C q and £ q of the Weyl algebra also satisfy the commutation 
relations of the s\.2 Lie algebra, namely 

(3.14) [n q ,C q } = -£ q , [£ q ,K q ]=2K q , [£ q ,C q ] = -2C q 

(which can also be checked by brute force calculation). This leads us to the following 
important result. 

COROLLARY 1. The q-ordering map O q : P(M. 2d ) — > W^d is an intertwining map for 
actions of the Lie algebra 5I2 realized by the triples {R, L, E} and {lZ q , C q , £ q }. 

4 The analysis of radial Weyl polynomials. We now apply the idea of R. Howe 
to use the action of 5I2 determined by {TZ q , £ q , £ q } to decompose the Weyl algebra W. 
Given any homogenous polynomial p G P l we can decompose it according to (2.4), so 
that p = Yl=l r 2k h k , h k e H l - 2k . It follows that 

[1/2] 

(4.1) o q ( P ) = Y,K°M- 

k=0 

By virtue of the intertwining property ( |3.10[ ) and (3.12) we obtain: 

Proposition 1. The image of the space H of harmonic polynomials in P under the 
q-ordering map O q is independent of q. For w € W we have 

d 

weO q {H) <=► ]T[a„[ a +, W ]] = 0. 
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We shall denote this space by H = O q (H) and shall refer to it as the space of Weyl 
algebra harmonics. We note that Geller in Q observed, with a very different method and 
for the particular cases q = 0, |, 1 only, this somewhat surprising fact that the space of 
Weyl algebra harmonics does not depend on the ordering. 



Getting back to (4.1), we see that two factors are essential for understanding the 



structure of the Weyl algebra: the space Tt C W of Weyl harmonics and the action of 



operators lZ q on W. In fact, the equality (4.1) leads to the following result 



Proposition 2. Recall the space of radial polynomials in P, J — C[r 2 }, and for 
Pi 1 " 2 ) = J2k=o c k r2k € J set p(TZ q ) = J2k=o c kK q . Then the map 

JxH3(p,h) — ► p{K q )O q {h) G W 
gives rise to an isomorphism of J ® H with W. 

In the sequel we shall investigate in detail the image J q — O q (J) — {p(lZ q )I \ p G J} 
of the space J of invariant polynomials under the ordering map. Although apparently 
depending on q through its construction, it will turn out to be the same for all values of 



q and consists of polynomials of a single element (the number operator (4.3)) of W — cf. 
Proposition || below. 

For each nonnegative integer k we set 

(4.2) r, k =r, k [q,d]=n k q I, 

so that {rjk} is a natural basis of J q . Let us write down the explicit f orm of a few first 
elements of this sequence. Clearly 770 = I, and after collecting terms in ( |3.11 ) and making 
use of the CCR ([O]) we arrive at 

d 

m = Y^ a t a i + d ^ -q)i = N + d ( l - 

=1 



Here we have denoted by 



(4.3) N = Y,aj aj 

3=1 

an element of the Weyl algebra, in physics called the number operator, which plays a 
rather special role in questions related to the structure of the Weyl algebra, see e.g. 
|^|, Its relevance for the present situation stems from the fact that the space J q 

consists of polynomials of N. Before proving this fact, we bring in for illustration the 
expression of 772 in terms of N. Setting for simplicity to = c?(l — q), one reads it as follows 

m {N) = (N + t /) 2 + (1 - 2q)(N + t a I) + qt I. 



Lemma 4. For each k G Z + there exists a unique polynomial LUk{t) = LUk(q,d;t) of 
precise degree k, LUk G C[t], such that 

r]k[q, d] = ujk(N). 

Proof. Clearly u>o(t) = 1 and the form of uJi(t) and w 2 (£) can be read of the above 
formulae. To prove the result in general we proceed by induction with respect to k. First 
note the general relations (the so called pull-through relations) 

(4.4) a l p{N)=p(N + l)a l , a+p(N) = p(N - l)a+ 
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valid for any polynomial p(N) of the number operator. It follows that 

(4.5) TZ qP (N) = q(l - q)(2N + d)p(N) + (1 - q) 2 (N + d)p(N + 1) + q 2 Np(N - 1). 

A comparison of highest order coefficients in terms at the right hand side shows that 
applying lZ q to a polynomial in N increases its degree by one. □ 
For completness we explicitly state the following 

Proposition 3. The space J q is independent of q and consists of polynomials in 
terms of the number operator N . Thus we set 

(4.6) J = J q = C[N]. 

To determine the form of polynomials Uk we shall investigate the difference and re- 
currence equations satisfied by them. We begin with the following observation. 

Lemma 5. The s pace J is invariant under the action of the triple (7Z q , C q , £ q ) and 
in fact rjk given by (4.2), with k ranging over all nonnegative integers, form a weight 
basis for the sfe action on J defined by it. Setting for convenience rj^i — this action is 
expressed by the well-known formulae 

TL q r\ k = r]k+i, 

£ q r)k = k(k + d- l)?^-!, 

£ q n k = (2k + d)rj k . 



This immediately follows by comparing the definition (4.2) and the corresponding action 
of {R, L, E} on the space of radial polynomials. 

By virtue of Lemma ^ we can transfer the action of {1Z qi C q , £ q } to the space of 
polynomials C[t]. If X is any from the triple of operators {7Z q , C qi £ q } we define the map 
X : C[t] i-> C[t] by setting 

(4.7) {Xp){N) = Xp{N). 

By the proof of Lemma ^, the right hand side is a uniquely determined polynomial of the 
number operator N. 

Employing the classical notation for the forward and backward difference operators, 
i.e. setting 

Aw(t) = u>(t + 1) - Lu(t) Vw(t) = - u>(t - 1), 
we can state the following. 

Lemma 6. The action of the triple {TZ q , C q , £ q } on J induces the following action 
on C[t}. 

R q uj{t) = {l-q) 2 (t + d)Au>(t)~q 2 Woj(t) + (t + d(l-q))uj(t) 

Z q oj(t) = (t + d)Au(t)-tVuj(t), 

Equ(t) = 2(1 - q)(t + d)Aw(t) + 2qt\7uj(t) + duo(t). 

The proof is a routine calculation based on Lemma || and (4.4). 
A consequence of the relation 

£ q r 1k = E q uj k (N) = ku>(N) 
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is a difference equation satisfied by lo^ which reads 

(4.8) qtVcu k (t) + (1 - q)(t + d)Aw fc (t) = ku k (t). 

In order to fix a multiplicative constant for io k which is not determined by the equation 



(4.8) , we look for an additional relation. This is supplied by taking a suitable combination 
of operators appearing in Lemma|. One notes that the operator R q — q(l— q)L q +(q— \)E q 
acts on polynomials of t as multiplication by the factor (t + |)- Thus we see that the 
polynomials lou satisfy the recurrence formula: 

(4.9) Lo k+1 {t) - q{l - q)k(k + d- l)u k -i(t) - [t + (1 - q)d - (2q - l)k]u k (t) = 0. 



Setting in addition w_x = 0, the formula (4.9) makes sense for all k £ Z + (= {0, 1,2,.. .}) 



and with the initial condition ujq — 1 determines the sequence cut uniquely. Now we can 
formulate our first main result. 

Theorem 1. Let 2 F\(a, 6, c; x) be the Gauss hypergeometric function. The polynomials 
uju defined in the Lemma || are given by the formula 

Wk(t) = (rf)fe(l - q) k 2F 1 (-t, -k, d; )> 

where (d)k = d(d + l)(d + 2) . . . (d + k — 1) denotes the shifted factorial. 

Proof. Let us first recall the well known Gauss recurrence formulae for the hypergeo- 
metric function (cf. 

(2b — c — bx + ax) 2 Fi(a, 6, c; x) + (c— b) 2 F 1 (a, 6—1, c; x) + b(x— l) 2 Fi(a, 6+1, c; x) = 0. 
After the substitution of parameters 

a = —t, 6 = —k, c = d, 



l-q" 



it coincides with the recurrence equation (4.9). Since by definition ujq — 1 and similarly 
2 Fi(—t, 0, d; jz-q-) = 1 we obtain the required result. □ 

Remarks. 

1) For q = 1 the expression for uj^ given in Theorem 1 should be understood as a limit 

fc - . • • - 1 



w fc (t) = (d) k lim (1 - qy 2 F 1 {-t, -k, d- 



q^l- ' l-q 

2) Other contiguous relations for 2 Fi give other formulae fulfilled by lo^. For example 
from the equation 

(2a — c — ax + bx) 2 Fi(a, 6, c; x) + (c — a) 2 Fi(a — 1, 6, c; x) +a(x— l) 2 Fi (a+ 1, 6, c; x) = 



we obtain the difference equation ( J4.8|) 



With a suitable change of variable it is possible to obtain an explicit expression for 
the generating function for this family of polynomials. Recall that given a sequence of 
polynomials (fffe(A))^L , the function 

oo 

G(\,s) = Y.9kW sk 
is said to be a generating function for (<?fc(A))^ . 
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Now make the substitution t + to = — —A, where a = (q(l — q))^ 1 ^ 2 and to = (1 — q)d 
and suitably renormalize the sequence of polynomials by setting 

(4.10) g k (X) = g k {q,d; A) = —^-uj k {t), where A = ia(t + t Q ). 



Setting in addition sq = ia(q — i), we can rewrite the recurrence equation (LE) in the 
form 

(4.11) (k + 2)g k+2 {\) + (k + d)g k {\) + [-A + 2s (fc + l)] 5fe+ i(A) = 0. 



Multiplying both sides of (4.11) by s k+1 and taking the sum over all nonnegative k, we 
obtain an equation satisfied by the generating function for (<7fc(A))J*L : 

BC 

(l + 2s s + s 2 ) — (X,s) = (-ds + X)G(X,s). 

OS 

Performing integration over s, we arrive at the following result. 

Theorem 2. Given q e [0, 1] and d e N the function 

exp [-(A + ds ) arctan (-(s + s ))) 

4.12 G(q, d; A, s) = P[aK °> J°} °ML 

[{s + s o y + ^]2 

has the following expansion 

oo 

G(q,d;X,s)^Y,9k(X)s k . 

k=0 

Now we we discuss in more detail the classical case q — \ of the Weyl ordering. By a 
straightforward inductive argument we obtain 

Lemma 7. The polynowials w k (h,d;t) are even for even k, and odd for odd k. 



As the consequence, after the substitution ( 4.10 ), the polynomials gfc(A) are real valued 
as functions of A = ia(t + to). Moreover, the previous results simplify to the following 
form. 

COROLLARY 2. The sequence (g k ) satisfies the recurrence equation 

(4.13) {k + 2)g k+2 (X) = Xg k+1 (X) - (k + d)g k {X), 
and its generating function G(d] A, s) = d; A, 5) is given by 

„A arctan s 

(4.14) G(d;X,s) = —_—-, A; = 0,1,... 

(vs + 1) 

As we mentioned before, a special case (with d = 1 and q = |) of these results was 
obtained previously by C. M. Bender et al. in M, and more systematically by T. H. Koorn- 
winder in ]To{ . In particular the connection of g k with the continuous Hahn polynomials 
is observed there, and in the first mentioned paper also the formula for the generating 
function identical with the case of d = 1 of our formula (4.14) was given. However, neither 
the case of q = \ and d > 1 nor other cases of our results were mentioned there. 

The next result identifies the polynomials g k (X) with known classes of orthogonal 
polynomials. 
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Recall (see jlO]]) that the continuous Hahn polynomials p k (x; a, b, c, d) are defined by 

Pk(x; a, b, c, d) = i k - — — — — ~ — — 3-p2(— fc 5 fc + a + 6 + c + c? — 1, a + ix; a + c, a + d: 1) 

k\ 

where 3F2 is generalized hypergeometric function (see 0), while the Meixner-Pollaczek 
polynomials P^?\x; <j>) are given by 

pW(a» 0) = e ^ 2 F! (-n, a + tar, 2a; 1 - e -34 *). 

The proof of the next result is a matter of standard calculation, when one takes into 
account the following formula connecting hypergeometric functions, quoted in |To| . 

2F 1 (-n, 2a + 2ix,Aa; 2) = 3 F 2 {-n, n + 4a,a + ix; 2a, 2a + -; 1). 



Theorem 3. For an arbitrary d £ N and for o—\ the polynomials g'fc(A) = <7fc(|, d; A) 
given by (4.10) can 6e identified as follows 

(d) fc /Add 1 d d 1 



(Dfcd + Dfe V4'4'4 2' 4' 4 2 

and a/so 

Moreover, they are orthogonal on (—00, 00) with respect to the measure p(X)d\, where 
p(X) = \T(l + *f)| 2 . 

Let us point out that the polynomials LUk(t) do not share certain properties of g k (X). In 
particular by rewriting the recurrence equation ( |4.9| ) in the form 

<*>fc+i(i) = (^fc* + B k )uj k {t) - C k u k -i(t) 

we see that 

A^iAfcCfc = -g(l - q)k(k + d- 1) < 0. 
Thus the necessary condition for the existence of a positive measure, with respect to 
which ujk would be orthogonal, is violated. 

Also in the general case of nonsymmetric ordering corresponding to q 7 ^ ^ the general 



picture is not so clear and satisfactory as in the symmetric case. From (4.8) one could 



expect that the set {u k } can be described as a subset of two known families of polynomials 
— Krawtchouk (K n ) and Meixner (M n ). In fact by comparing parameters we see that 

2 F 1 (-t, -k, d; ^—) = K k (t, 1 - q, -d) = M k (t, d, -l—l). 

However, there is again a problem with orthogonality, since the parameters occuring in 
these formulae fall outside the range required by the orthogonality condition (see JIJ and 
0). 
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